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Using archival observations recorded over a 5+ year time- 
frame with the NASA Galaxy Evolution Explorer (GALEX) 
satellite, we present a study of the ultraviolet (UV) variability 
of 4360 quasars of redshifts up to 2 = 2.5 that have optical coun- 
terparts in the Sloan Digital Sky Survey DR5 spectroscopic cat- 
alog of Schneider et al. (2007). The observed changes in both 
the fai- UV (FUV: 1350 - 1785A) and near UV (NUV: 1770 - 
2830A) AB magnitudes as a function of time may help differ- 
entiate between models of the emission mechanisms thought to 
operate in these active galaxies. A list of NUV and FUV vari- 
able quasars was derived from the UV light-curves of sources 
with 5 or more observational visits by GALEX that spanned a 
time-frame > 3 months. By measuring the error in the derived 
mean UV magnitude from the series of GALEX observations 
for each source, quasars whose UV variability was greater than 
' ^' the 3-cr variance from the mean observed value were deemed to 
I be (intrinsically) UV variable. This conservative selection crite- 
_ rion (which was applied to both FUV and NUV observations) 
. resulted in identifying 550 NUV and 371 FUV quasars as being 
^ • statistically significant UV variable objects. 
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1. Introduction 

The flux emitted from quasars and other Active Galactic Nuclei 
(AGN) such as Seyferts, BL Lacs and Blazars has been found 
to vary aperiodically on time-scales from hours to decades and 
from X-ray to radio wavelengths. In the case of radio-loud and 
radio-quiet quasars optical flux changes over timescales as short 
as minutes and with amplitudes ranging from a few hundredths 
to a few tenths of a magnitude (often refereed to as 'optical 
micro variability'), have been widely reported (Miller, Carini & 
Goodrich 1 19891 Stalin et al. '2005V The mechanism(s) respon- 
sible for this short-term microvariability is still debated, with 
much of the visible spectral variability being explained by ther- 
mal processes associated with a black-hole accretion disk and/or 
non-thermal processes associated with the relativistic beaming 
of a jet (Ramirez et al.'2010', Pereyra et al.'2006'). In fact Sesar et 
al. ( 2007 ) have found that nearly all quasars vary at some level 
when observed over long enough times scales, with ~ 90% of 
quasars found in the Sloan Digital Sky Survey (SDSS) varying 
by > 0.03 mag over a six year time scale at visible wavelengths. 
Although the actual cause of (long-term) variability in these ob- 
jects is also still a matter of intense debate, variability studies 
of these sources are important since (in principal) they can set 
limits on the size of the emitting regions and can also constrain 
models of the emission mechanisms thought to operate in these 
active galaxies. 

Numerous variability studies have been aimed at establish- 
ing correlations between quasar visible variability and important 
physical parameters such as source luminosity, redshift, time-lag 
(i.e. the time period between two sets of observations) and radio 
loud ness (H elfand et aL lIOOTI V anden Berk et al. l2004l de Vries 
et al. 12005 ', Ramirez et al. |2009l l. The establishment of such cor- 
relations is important since there is increasing evidence that the 
source of quasar variability is linked to the accretion disk sur- 
rounding the central supermassive black hole engine (Peterson 
et al.E M Pe reyra et aL l2006l Wold et aL l2007b . For example, 
Bachev ('2009) has studied the observed time lags between the 
continuum emission from several optical bands for a large sam- 
ple of quasars. He demonstrates that the lags are broadly con- 
sistent with an emission model in which the optical variations 
are largely due to the reprocessing of the central X-ray radia- 
tion in a surrounding thin accretion disk. However, both Kelly et 
al. (2009) and MacLeod et al. (2010), using a 'damped random 
walk model', have suggested that these optical emission varia- 
tions result from thermal fluctuations in the accretion disk that 



are powered by an underlying stochastic process such as a tur- 
bulent magnetic field. Clearly there is still much debate as to 
the interpretation of these optical variability studies (Kozlowski 
et al. 2010), and other emission mechanisms such as bursts of 
supernovae events close to the nucleus and extrinsic variations 
due to line-of-sight microlensing events cannot (as yet) be to- 
tally discounted (Ulrich et al. |1997| ). It is mostly agreed that the 
level of quasar variability increases with decreasing rest wave- 
length, which is related to the observation that their spectra be- 
come bluer in their bright phases. Several authors have claimed 
a variability-redshift correlation in the visible, although the ac- 
tual sign of the correlation in the studies is debated (Giveon et al. 
1999, Trevese and Vagnetti 2002). For example, Vanden Berk et 
al. (2004) have claimed evidence for redshift evolution in quasar 
variability, with greater variability occurring at higher redshifts. 
However, such a trend is only revealed when an allowance is 
made for rest wavelength, absolute magnitude and time lag for 
each redshift interval. Such corrections to the observed optical 
variability data are crucial, since spectral features may move in 
and out of the photometric pass-bands with changes in the target 
redshift. 

Quasars and AGN have also been well studied at ultraviolet 
(UV) wavelengths, such that typically their far UV spectra ap- 
pear as several strong emission lines (formed in gas located fur- 
ther away from a central black hole and excited by continuum 
radiation) superposed on an underlying accretion-driven contin- 
uum. In contrast, their near UV spectra contain only a UV con- 
tinuum with a 'bump' at 2800A superposed with the Mgll emis- 
sion lines (Vanden Berk et al. 2001). The UV variability of these 
objects has been widely reported, particularly with reference to 
reverberation mapping observations that in principal allow mea- 
surement of the size of the broad emission-line region and can 
infer the mass of the alleged central black hole (Peterson et al. 
1993, Kaspi et al. 20001 l2007l l. UV photons are generally be- 
lieved to be radiated from the quasar supermassive black hole 
accretion disk when its gravitational binding energy becomes 
released, while X-rays are produced in the disk corona. A sta- 
tistical survey of the UV variability of 93 quasar/ AGN, as ob- 
served with the I.U.E. satellite, was carried out by Paltani and 
Courvoisier (.1994 ) who found no strong evidence for correla- 
tions between variability and parameters such as redshift, lumi- 
nosity or UV spectra-index. However, changes and instabilities 
in the accretion disk must surely be linked to their observed UV 
variability, but recognizing and understanding the emission sig- 
natures from such features still remains problematic. 
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To proceed further with studies of UV variabiHty we require 
access to far larger data sets than previously possible, and with 
the launch of the NASA Galaxy Evolution Explorer (GALEX) 
satellite (Martin et al. 2005) we now have a significantly large 
database of far UV (FUV: 1350 - 1785 A) and near UV (NUV: 
1770 - 2830A) observations of quasars and AGNs that cover a 
timespan of over 5 years. Trammell et al. (2007 ) have presented 
a preliminary analysis of the basic UV properties of 6371 SDSS 
Data Release 3 (DR3) selected quasars observed during the first 
few years of the GALEX mission. One major conclusion from 
this study was that > 80% of the DR3 listed quasars had a cor- 
responding GALEX NUV detection, implying that many more 
quasars should be detected during the remaining years of opera- 
tion of the GALEX satelUte. Bianchi et al. ( I2009I I have shown 
that GALEX is also ideal for detecting quasars with blue (FUV- 
NUV) colors that are difficult to identify from visible observa- 
tions alone. Since the publication of these initial studies, newer 
and far larger versions of both the SDSS and GALEX catalogs 
are now in-hand. In particular, the catalog of quasars contained 
in the SDSS Data Release 5 (DR5) published by Schneider 
et al.( 120071 ) lists some 77,000-1- objects (with associated red- 
shift and luminosity values), the majority of which should have 
UV counterparts observable with GALEX. Furthermore, since 
many quasars are known to be variable, this potentially enables 
us to study the UV variability of a significantly large number of 
these sources. 

The largest group of known variable sources listed in the 
GALEX UV Variability-2 catalog are active galaxies (Wheatley 
et al. 120081 ). This sample of sources was derived solely through 
selection due to their observed UV variability, without prior 
knowledge that such sources were in fact AGN. These variable 
sources possessed redshifts as high as 2 2.5 and flux changes 
as large as ~ 2.3 AB mag. were found in both of the GALEX 
near and far UV bands for AGN that were repeatedly observed 
over a 3.5 year time period. Based on the apparent quality of 
these data, we decided to search the entire GALEX archival 
database (which now spans 5+ years of GALEX observations) 
for spatial matches to quasars previously identified and listed in 
the SDSS DR5 catalog and then, for those sources that had mul- 
tiple visit measurements, determine which ones exhibited mea- 
surable UV variability. This Paper discusses which GALEX 
sources were identified with SDSS quasar counterparts, which 
ones showed long-term UV variability (as measured over time- 
frames > 3 months) and which, if any, correlations exist between 
this UV variability and other astrophysically interesting source 
parameters. 

2. Data Analysis 

In order to produce a list of quasars that was suitable for an ul- 
traviolet variability study we first searched the GALEX GR5 
data archival release (at the Multi-Mission Archive at the Space 
Telescope Institute, MAST) for sources that spatially matched 
optically identified quasars listed in the SDSS DR5 spectro- 
scopic catalog of Schneider et al. (120071 ). A spatial match was 
deemed positive if the central position of the matching GALEX 
source was within ± 8 arc sec of the listed SDSS position, 
this number being slightly larger than the positional error of the 
GALEX observations for a point source detection. This search 
resulted in 77, 249 quasars that were common to both the SDSS 
and GALEX ultraviolet source catalogs. 

In order to produce a list of quasars whose UV variabil- 
ity could be accurately assessed as a function of time, we then 
down-selected the large list of UV-selected quasars using the 



following criteria. Firstly, we only selected quasars with either 
(i) GALEX observations in the All-sky Imaging Survey (AIS) 
that were brighter than NUVmag = 19.0 and with an exposure 
time of at least 100s, or (ii) those quasars with GALEX Deep 
Imaging Survey (DIS) observations of 1500s duration or more 
(see Morrissey et al. l2007l for a detailed discussion of the obser- 
vational modes of the GALEX satellite). A further restriction 
was applied such that only those quasars that had 5 or more ob- 
servational visits by GALEX (and were detected in at least one 
of those visits) were selected. We note that an observational visit 
by GALEX does not necessarily result in a detection, since an 
intrinsically variable object may fall below the detection limit 
throughout a particular visit. We were guided in these choices by 
Figure 4 of Trammel et al. (2007 ) which shows the typical rms 
magnitude errors for an AIS and DIS observation with GALEX. 
Those results demonstrate that the measurement error for a typi- 
cal AIS observation of a relatively UV bright (NUV„ia<, = 19.0) 
object is approximately equal to the measurement error for a 
1550s DIS observation of a faint (NUVmag > 21.0) source. In 
Figure 1 we show a plot of the number of these selected quasars 
as a function of the number of observations in the GALEX NUV 
and FUV channels. It can be seen that there is a significant num- 
ber of sources that have been visited between 5 and 20 times, 
with several sources being repeatedly observed > 100 times. On 
application of the selection criteria of both exposure time and >5 
visits, this resulted in a list of 4360 quasars that were detected in 
at least one observational visit by GALEX. 

In Figure 2 we show a plot of these selected quasars as a 
function of the timespan covering all of their 4360 NUV and 
2709 FUV GALEX observational visits (NB: we note that the 
sample of FUV selected quasars is smaller than the number of 
NUV selected ones due to operational constraints on the FUV 
detector that often precluded simultaneous FUV and NUV ob- 
servations of a source). This Figure demonstrates that our quasar 
sample is not biased towards any particular timeframe of vari- 
ability, with the majority of selected sources having been repeat- 
edly sampled over a 1 to 5 year timespan. We note that although 
the time between many observational visits to a particular quasar 
may span several years, the target may be repeatedly observed on 
time-scales < 1 day during a series of observational visits with 
GALEX. 

All the data shown in Figures 1 and 2 clearly indicate that our 
sample of quasars is both large and contains significant numbers 
of objects that have multiple (> 5) observational visits spread 
over a significantly long timespan. In order to ensure that quasars 
with potentially short-term variations of a few days had > 5 ob- 
servations, we have further restricted our target selection to those 
ones whose visits spanned a time-frame > 3 months. This se- 
lection criterion resulted in a final list of 3521 potentially UV 
variable quasars, of which 1942 had 10 or more observational 
visits. To ensure that all of the UV-selected sources were not in 
any preferential direction on the sky, in Figure 3 we show the 
galactic locations of our SDSS selected sample that have been 
observed 5 or more times. Since our selection of targets was 
derived from the (northern hemisphere) SDSS survey catalog, 
clearly only quasars with declinations > -12° are available for 
present study. Although there are large areas of the sky where 
we have no UV-selected quasars. Figure 3 shows no obvious 
sight-line preference apart from most sources being observed 
at high galactic latitude away from the effects of absorption of 
ultraviolet flux by gas and dust in the galactic plane. We then 
searched this list of 3521 sources for the ones which possessed 
statistically significant UV variability between each of their ob- 
servational visits. Although the quasar nature of these sources at 
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visible wavelengths was already known from the SDSS DR5 cat- 
alog, we decided that statistical confirmation of their associated 
UV variability was also required for this study. Due to the limited 
sensitivity of the GALEX observations this naturally restricted 
our study to the more UV variable quasars in this sample. 

2. 1. UV Variability Selection 

In order to determine which of the 3521 selected quasars were 
variable at least once during their series of observations, we re- 
quire knowledge of the size of the error associated with the es- 
timation of each of the NUV and FUV magnitudes measured 
for each visit to a particular source. This error in this measure- 
ment is affected by the length of the observation by GALEX, 
the position of the source on the UV detector and the associated 
varying sky background signal which are both dependent on the 
(changing) orbital viewing characteristics of the GALEX satel- 
lite. Thus, in order to progress further in our study we clearly 
require a more consistent estimation of the measurement error 
for the UV magnitudes derived from each observational visit by 
the GALEX satellite. Fortunately, Gezari et al. (|2008Jl have es- 
tablished a robust numerical method of assessing the statistical 
significance of the variability of a series of GALEX UV ob- 
servations of the nuclei of normal quiescent galaxies, which is 
directly applicable to our present QSO data set. The method in- 
volves deriving the quantity, fT(<m>), which is defined as the 
standard deviation (a) from the weighted mean (<m>) of a se- 
ries of observed UV magnitudes (m). For a series of observa- 
tions of the same object, the quantity of a{<m>) will vary such 
that above some statistical level any observed change will be due 
to intrinsic source variability. Gezari et al. (2008) estimated the 
Poisson errors for a series of UV magnitude determinations us- 
ing GALEX (see equation (2) in that paper). They found that 
the measured standard deviation from a mean UV source magni- 
tude was dominated by systematic errors that were not accounted 
purely by the photometric error We note that the major contri- 
bution to the error in the estimation of a UV magnitude using 
GALEX is thought to be due to the (often varying) sky back- 
ground level. Thus to be conservative, Gezari et al. (2008) used 
a cut of 5-a above the mean magnitude level to determine intrin- 
sic source variability. 

Values of cr(<m>) for our quasar sample were also derived 
from equation (2) of Gezari et al. (I2008I I for each set of obser- 
vational visits. The values of NUV and FUV AB magnitudes 
recorded for each visit to a selected quasar, together with their 
measurement errors, were retrieved from the GALEX GR5 data 
archival release at the MAST. We used circular aperture magni- 
tudes of a 6 arc sec radius, as listed in the MAST catalog, which 
were obtained using the GALEX pipe-line generated SExtractor 
magnitude estimator (Bertin & Arnouts [19961 1 . In Figure 4 we 
show a plot of a(<m>) versus <m> for the multiple NUV ob- 
servations of 3521 of the selected quasars. As one might expect, 
the value of <j(<m>), which is essentially a measure of the er- 
ror in the derived mean AB magnitude of a quasar, increases with 
the faintness of the object under study. All of these 3521 quasars 
were detected as being optically variable in the SDSS DR5 cat- 
alog, and although they also all may be UV variable (though this 
has never been actually proven in any large study), we have de- 
cided to take a conservative approach in selecting only the more 
UV variable objects for further study in this paper To this end 
we have decided to use a threshold cut using a standard deviation 
of 3-cr above the level of cr(<m>) for a given mean magnitude 
(<m>) to reveal significant UV variability. On application of 
this threshold cut we find that the smallest detectable magnitude 



variation is 0.1 1 mag, which was determined for a 16.5 AB mag 
quasar of redshift z = 0.4. 

The 3-cr threshold cut is shown superposed on both Figures 4 
and 5 respectively for the NUV and FUV observed magnitudes. 
Application of this statistical significance selection criterion re- 
sulted in the identification of 550 quasars as being NUV variable 
and 371 quasars as being FUV variable, with the vast majority 
of the FUV variable objects being common with those identified 
as variable in the NUV. For simplicity, from now on we term 
this sub-set of selected variable objects as UVQs (i.e. ultravio- 
let variable quasar stellar objects). The typical levels of observed 
UV variability for each of these UVQs, can be best viewed when 
shown as NUV and FUV light-curves (i.e. the AB magnitude 
for each visit versus observation date). In Figure 6 we show 
the FUV and NUV light-curves for a relatively bright (SDSS 
J163915.80H-412833.6) and a faint (SDSS J221629.33H-002340) 
quasar, that were both repeatedly observed over a 4+ year time- 
frame with GALEX. It should be noted that faint and variable 
objects can fall below the GALEX threshold of detection during 
a series of observations, and hence the derived mean magnitude 
(which is computed from only actual detections within a series 
of visits) will be of a higher value than had the object been de- 
tected during all visits. Thus, the computed change in magnitude, 
which is a measure of the object's variability, will be underesti- 
mated. Finally we note that we did not discover any UVQs that 
were highly variable in the FUV but not variable in the NUV 

2.2. The Structure Function 

The mean variability of UVQs is often described by a Structure 
Function, which is the correlation between two points in a set of 
measurements of an object's light curve (Simonetti et al. |19851 l. 
Thus, in its simplest form, it provides a statistical measure of 
variability as a function of the time-lag between many sets of 
observations. However we note, as several authors have pointed 
out, there are limitations in its applicability to derive timescales 
and amplitudes of quasar variability (Vio et al. I1992I I. In addi- 
tion, quasar variability is most probably a complex function of 
several (possibly interdependent) physical parameters such as lu- 
minosity, rest-frame wavelength, redshift, black hole mass and 
time-lag. Helfand et al. (2001) have pointed out that 'the com- 
plex interdependence of the observed variability (of quasars) on 
wavelength, redshift, bolometric luminosity, radio loudness and 
possibly other parameters makes it difficult to draw definitive 
conclusions regarding the physics of the central engine or quasar 
evolution'. Thus, 'caveat emptor' should be the watchwords for 
any discussion of quasar variability dependencies, and thus we 
employ a method by which these dependencies of variability can 
be plausibly disentangled. 

For our present study, in which most of the identified UVQs 
each have > 5 measurements, we start by using the standard 
formulation of the Structure Function (SF) given by di Clemente 
et al. (lT996l l: 

SF= VV2h<(| Am(Ar) IP - (a^) 

where Am(Ar) is the difference in UV magnitude be- 
tween any two observations of a UVQ separated by Ar in the 
quasar rest time frame, and (t„^ is the square of the uncertainty 
in that difference (being equal to the sum of the two individual 
observations' errors in quadrature). The final SF for all of our 
UVQ observations is just the rms of the summation for all 
possible combinations of magnitude measurements for a given 
value of time-lag. 



Welsh, Wheatley & Neil: Quasar UV variability 



5 



In Figure 7 we show log-log plots of the SF (i.e. the mea- 
sure of UV variability) versus the rest frame time lag (in days) 
for both the NUV and FUV observations of each of the statisti- 
cally significant variable UVQs selected from Figures 4 and 5. 
Following the work of Vanden Berk et al. (i2004j we have plot- 
ted the data averages in 10 time-lag bins, chosen with equal in- 
tervals in the logarithmic rest-frame time lag. Note that the data 
as presented in Figure 7 make no account for any other plausi- 
ble variability dependencies such as luminosity, redshift or radio 
loudness. In this form, the variability function is often termed 
an 'ensemble' Structure Function. EiTor bars were determined 
by propagating the rms errors derived for the average magnitude 
difference and measurement uncertainty in quadrature for each 
observation. 

The data shown in Figure 7 reveal two important facets of the 
UV variability of quasars. Firstly, the amplitudes of variability 
in both the NUV and FUV are far larger than those observed at 
visible wavelengths (which are typically only ~ 0.03 mag rms), 
and secondly, the levels of FUV variability are greater than those 
observed in the NUV. 

The NUV SF shown in Figure 7 increases linearly (with a 
slope of 0.439) from a rest frame time-lag of 1 to ^100 days, 
with a definite turn-over and subsequent flattening for time-lags 
^ 300 days. The actual best-fit to all of the NUV data is a 
quadratic in log-log form. A similar correlation of variability 
with time-lag has also been found in many other studies, and 
the flattening of the curve at long time-lags has been reported 
by both Vanden Berk et al.( |2054] | and C ristie et al. (fT997] l. 
However we note that de Vries et al. (I2005I I. in an optical survey 
of long-term quasar variability up to timescales of ^ 40 years, 
has observed that this apparent flattening of the SF only lasts for 
a few years before continuing its linear increase with time-lag. 
In essence, they claim no significant turn-over of their SF, such 
that visible quasar variability increases monotonically with in- 
creasingly long time-lags. Although our sampling of UVQs is 
limited to time-scales of < 5 years, inspection of the behavior of 
the NUV SF curve suggests that the rate of increase in UV vari- 
ability as a function of time decreases when observed over time- 
lags greater than ^ 2 years. The FUV Structure Function, also 
shown in Figure 7, follows the same general trend as the NUV 
curve, except that we have observed a greater level of baseline 
variability of UVQs in the FUV. However, the change in FUV 
variability as a function of time-lag is less pronounced than that 
in the NUV, as indicated by the best-fit line of slope 0.292. The 
FUV data, like those recorded in the NUV, also show a flattening 
in UV variability for time-lags ^ 300 days. 

3. Variability Correlation Plots 

Previous studies of quasar variability in the rest-frame of the 
optical/UV regime have revealed several apparent physical cor- 
relations. For example, trends of variability with wavelength 
and luminosity (first discovered by di Clemente et al. 1 19961 and 
Hook et al. ,1994, respectively) have been confirmed by Vanden 
Berk et al. (I2004I I who have found that quasars are more vari- 
able at shorter rather than longer wavelengths and the ampli- 
tude of variability decreases with increasing luminosity (i.e. 
intrinsically brighter objects vary less). These results, and the 
SF method used to obtain them, have recently been verified 
for a sample of quasars observed over a 10 year period by 
MacLeod et al. (120101 1. Furthermore, radio-loud quasars appear 
more variable than their radio-quiet counterparts (Helfand et al. 
1200 11 1, and the apparent increase in variability for higher red- 
shift quasars is mainly due to the increase of variability with 



rest-frame frequency and to the increase of rest-frame frequency 
scanned at high redshift for a given frequency (di Clemente et 
al. 1996). The Structure Function derived in the previous sec- 
tion, which described UV variability as a function of rest-frame 
time-lag, took no account of any of the other possible inter- 
dependencies such as luminosity or redshift. To reveal such multi 
inter-dependenices would require 4-dimensional plots, and thus 
we need an alternative method to explore other possible correla- 
tions with the UV SF. Therefore, we closely follow the work of 
Vanden Berk et al. ('2004') who have described a method that at- 
tempts to disentangle these possible physical interdependencies. 
Since each of the 10 time-lag bins of the UVQs plotted in Figure 
7 contains objects with a wide-range of other (differing) phys- 
ical parameters, we start by subdividing the targets contained 
within each time-lag bin into separate (and assumed to be de- 
generate) physical parameters. This then allows us to obtain a 
set of (quasi-independent) variability relations with respect to 
the SF for a particular time-lag, UVQ luminosity, redshift and 
radio loudness. In the following subsections we present these 
variability relationships, and then briefly discuss their possible 
astrophysical meaning in Section 4. 



3. /. Time Lag - UV Variability Dependence 



Firstly we explore the dependence of UVQ NUV and FUV vari- 
ability on the rest frame time-lag, in a form that is (quasi) in- 
dependent of the parameters of luminosity and radio loudness. 
We divide all of the UVQs in our UV variability sample, whose 
redshift {z) values span the range of to 2.5, into three red- 
shift bins of size < 2 < 0.47, 0.47 < z < 1.33 and z > 1.33. 
The physical relevance of these 3 redshift bins with respect to 
observations through the two GALEX pass-bands will be dis- 
cussed in Section 4. In Figures 8 and 9 we show log-log plots 
of the respective NUV and FUV Structure Functions versus rest 
frame time-lag, with the data averages being plotted separately 
for each of these three redshift bins at any given rest frame time- 
lag value. Both of these plots, in agreement with the simpler en- 
semble SFs shown in Figure 7, show that the amplitudes of both 
the NUV and FUV variability of UVQs increase as a function 
of rest frame time-lag, irrespective of the value of quasar red- 
shift for time-lags < 200 days. For time-lags > 300 days there 
is a pronounced rollover in the NUV SF for all redshift values, 
which is also observed in the lower S/N FUV data. To further 
demonstrate this roll-over effect we have also included separate 
linear best-fits for each of the 3 redshift bin data sets in Figures 
8 and 9. These straight-line fits show that rate of increase of the 
NUV SF up to a time-lag of ^ 200 days is greatest for sources 
of redshift 0.47 < z < 1.33, an effect that is also miiTored in the 
FUV data set of Figure 9. For time-lags > 300 days it is clear 
that the same linear fits are inappropriate for all of the NUV and 
FUV data, such that a far slower increase in the rate of the NUV 
and FUV SFs occurs. 

Figures 8 and 9 both show that for time-lags < 20 days there 
is a definite trend for the higher- 2 UVQs to be more NUV and 
FUV variable than their low-z counterparts. This trend has also 
been seen for visible quasars by de Vries et al. (2005). However, 
for time-lags > 100 days this trend mostly disappears such that 
both the FUV and NUV SF variability amplitudes are similar for 
all values of UVQ redshift when observed over long time-lags. 
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3.2. Luminosity - UV Variability Dependence 

Unlike redshift, absolute luminosity is an intrinsic physical prop- 
erty of a quasar and thus an investigation of UV variability (i.e. 
SF) as a function of this parameter should be less affected by 
other physical interdependencies. We make use of the values of 
luminosity given in Schneider et al. ( 12007i for each of our UV 
variable sources, and note that these values were calculated us- 
ing extinction corrected SDSS (absolute) i magnitudes and with 
the assumption that the quasar spectral energy distribution in the 
ultraviolet-optical can be represented by a power law. 

We start by separating the possible dependency effects of rest 
frame time-lag by dividing the data into four time-lag bins of 1 - 
10 days, 1 1 - 100 days, 101 - 500 days and > 500 days. In Figure 
10 we show log-log plots of both the NUV and FUV variability 
amplitudes (i.e. the SF's) as a function of (averaged) absolute 
magnitude in bins of 2 mag width for each of the aforemen- 
tioned time-lag bins. The best-fits to these data show that there 
is a slight decrease in both the NUV and FUV SFs as a func- 
tion of increasing luminosity for time-lags > 100 days. Similar 
(visible) luminosity plots have also been constructed by Vanden 
Berk et al. (I200'4D for a sample of 25,000 quasars in the SDSS 
survey, in which they similarly found that intrinsically brighter 
objects vary less than intrinsically fainter ones. However, for 
shorter time-lags we find the opposite behavior in which the 
NUV and FUV SFs increase as a function of increasing lumi- 
nosity (we note the apparent exception of the NUV SF for time- 
lags < 10 days, but this may be due to small number statistics 
with associated larger eiTors). The NUV and FUV data shown 
in Figure 10 also show that variability amplitudes are greater at 
longer time-lags than at shorter time-lags for all values of UVQ 
luminosity. This behavior is also revealed in the SFs of Figures 
7-9, which were derived independently of luminosity. This dif- 
ference in variability for both short and long time-lag values is 
discussed further in Section 4. 1 with respect to the work of Kelly 
et al. ( I2009l l. 

We note that there are no low-luminosity UVQs detected at 
high redshift in our sample, such that the low luminosity bins in 
both plots of Figure 10 are populated mostly by low-2 objects. 
Therefore, since the data in Figure 10 may still contain some 
hidden dependency on quasar redshift, in Figure 1 1 we re-plot 
the NUV and FUV variability data but now subdivided into bins 
of redshift. For convenience of plotting we have selected repre- 
sentative redshift bins of 0.5 < z < 0.895 and z > 2.028 for the 
NUV data, and bins of 0.5 < z < 0.895 and 1.4 < 2 < 2.03 for 
the FUV data, all plotted for the same two representative time- 
lag bins as used in Figure 10. For the NUV data in Figure 11 
the 4 curves of SF versus luminosity do not precisely follow the 
same trends as shown in Figure 10, although there is an overall 
general trend of decreasing SF value with increasing luminosity. 
This general trend of decreasing variability amplitude is also ob- 
served in the FUV data, except for the case of high z and small 
time-lags where an increase of SF variability with luminosity 
is observed. Since the general trends seen in both Figures 10 
and 1 1 are similar, we can deduce that redshift does not seem to 
have a profound effect on the UV variability- luminosity nature 
of these objects, and also that the shape of the UV SFs is only 
weakly dependent on quasar luminosity. 

3.3. Radio Loudness - Variability Dependence 

In Figures 12 and 13 we plot the SFs for the NUV and FUV 
variability of our quasar sample against their rest frame time-lag 
for both radio loud and radio faint UVQs. Based on their peak 



20cm radio flux density as measured in the VLA FIRST survey 
and using the absolute SDSS i magnitudes and distances listed 
in Schneider et al. (2007) to derive a measure of their optical 
luminosity, we compute values of the ratio of radio to optical 
luminosity in the rest frame (R*). Using equations (1) and (2) 
of Sramek & Weedman (1 19801 1 to derive R*, we use the widely 
accepted definition of R* < 1 for radio quiet sources and R* > 
1 for radio loud sources. Since many of our variable quasars do 
not have listed values for their radio fluxes, our NUV sample is 
restricted to 484 radio loud and 66 radio quiet sources and our 
FUV sample contains 327 radio loud and 44 radio quiet sources. 
Hence, the reader is reminded that we are dealing with small 
number statistics for the radio quiet source data in both Figures 
12 and 13. 

Firstly, inspection of both Figures clearly shows that the 
overall level of SF variability is greater in the FUV than the NUV 
for both radio loud and radio quiet UVQs. This is merely a re- 
statement of the behavior of the entire UVQ sample, as shown 
in Figure 7. However Figure 12 also indicates that radio loud 
sources are ^ a factor 2 more NUV variable as a function of 
time-lag than radio quiet sources, and the plot also shows a larger 
scatter in the shape of the NUV SF for radio loud than for radio 
quiet UVQs. However, there is only a marginal difference be- 
tween the best-fit slopes to these NUV SF data for both types 
of radio source. Therefore, our data would appear to support the 
notion that NUV variability is not strongly linked to the radio 
properties of UVQs, and as such, the physical mechanism re- 
sponsible for the NUV variability of radio detectable UVQs op- 
erates irrespective of their radio flux. 

The FUV data shown in Figure 13 have far fewer measure- 
ments and larger eiTors than the NUV data, but again the radio 
loud UVQs appear to be more variable as a function of time-lag 
(especially for periods > 100 days) than the radio quiet sources. 
The best-fit slope of 0.16 to the FUV radio loud data is signif- 
icantly different from that of the best-fit slope of 0.059 for the 
radio quiet data. This FUV variability difference is most pro- 
nounced for time-lags > 100 days. However, we note that this 
conclusion is based on only a small number of UVQs. We note 
that several authors have claimed that radio-loud quasars are 
(marginally) more variable than their radio-quiet counterparts 
(Eggers et al. .2000. Helfand et al. I2001I I. whereas Bauer et al. 
(I2OO9I 1 have found no evidence for any change in the optical 
variability amplitude of quasars with increasing radio loudness 
(luminosity). Several studies have found that optical microvari- 
ability, operating over time-scales of minutes to hours, occurs in 
quasars regardless of their radio properties (Ramirez et al. 20091). 
Our present NUV (and FUV) observations support the notion in 
which radio loud UVQs are only marginally more UV variable 
than radio quiet sources, irrespective of the time-frame of obser- 
vation. 

4. Summary Discussion 

Before briefly discussing any physical interpretation of the cor- 
relations shown in Figures 7 - 13 of Section 3, we first provide 
a short summary of the underlying UV properties of the UVQ 
sample. We also remind the reader that due to our conservative 
approach in the selection of UV variable quasars, these findings 
may only apply to the more UV variable of this class of object. In 
Figure 14 we show a typical quasar rest-frame spectrum (adapted 
from Vanden Berk et al.i2001jl, which shows the prominent Lya, 
CIV, CIII and Mgll emission lines superposed on a significant 
underlying UV continuum level, which for wavelengths blue- 
ward of Lya is heavily extinguished by the Lya forest of absorp- 
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tion lines. The entire UV energy distribution is often referred to 
as the 'big blue bump', and is thought to be dominated by the 
lO^'^^K thermal emission from a central black hole accretion 
disk. Based on the sensitivity and observational strategy of the 
GALEX satellite, the great majority of quasars that can be de- 
tected possess redshifts in the range < 2 < 2.5 (Bianchi et al. 
I2007I I. For quasars with redshifts in the range from 2 = 0.11 to 
0.47 the strong Lya emission line is present within the GALEX 
FUV passband, and for redshifts of 2 = 0.46 to 1.33 the Lya line 
is present within the NUV passband. In addition, the Lya forest 
of absorption lines (which considerably reduces the UV flux in a 
quasar spectrum) affects the entire GALEX passbands for red- 
shifts of z > 1.33. Therefore, any observed UV variability within 
either of the GALEX FUV or NUV passbands will be highly 
dependent on which of the emission/absorption/continuum fea- 
tures present in the quasar spectral energy distribution have been 
redshifted into, or out of, those two passbands and also which of 
these features are deemed to be most responsible for the physical 
origin of the observed UV variability. Therefore, our selection of 
3redshiftbins (0.11 < 2 < 0.47, 0.47 < 2 < 1.33 and 2 > 1.33) 
for the plotting of data in Figures 8 to 9 represent the appear- 
ance and disappearance of the strong Ly-a line in the FUV and 
NUV bands. Additionally, because of time dilation effects, high 
redshift objects are potentially observed for less time than low 
redshift ones. This will possibly affect any study of the time-lag 
on UVQ variability, which is why we have binned our time-lag 
data in the rest frame. 

Presently, explanations of the primary cause of quasar vari- 
ability can be broadly split into 3 classes: (i) gravitational mi- 
crolensing, (ii) discrete-event (Poissonian) processes and (iii) 
accretion disk instability models. Vanden Berk et al. (I2004I I 
have discussed these various processes (see the many references 
therein) and their visible data favor an accretion disk instability 
(or jet) cause of the observed optical variability. Since many of 
our present results (shown in Figures 7-13) generally agree with 
the visible observations of quasar variability by Vanden Berk et 
al. (120041) . we shall only address whether a model involving ac- 
cretion disk instabilities may explain our UV variability data. 
However, the reader should note that no theoretical model, irre- 
spective of its underlying physical mechanism, currently exists 
that can be directly compared with their (or our own UV) obser- 
vations. 

Long-term (7.5 yr) changes in the UV continuum level of or- 
der 25% - 125% have been reported for low luminosity quasars 
(Kaspi et al. 120001 1. which is about twice the continuum variabil- 
ity observed over a similar period for high luminosity (high-2) 
objects (Kaspi et al. 120071 ). Interestingly for this latter sample of 
(11) quasars, none showed Ly-a emission variations to a level 
of < 7% over the 6-year period of observations. A similar case 
of non-variability of the Ly-a line has also been reported for ob- 
servations of 3C 273 recorded over a 15yr timeframe by Ulrich 
et al (II993I 1. Several authors have claimed that quasar continuum 
variations are dominated by processes that are intrinsic to the ac- 
cretion disk (Peterson et al. .2004 ), and Trevese, Kron & Bunone 
(1200 11 1 have suggested that the observed spectral changes are 
consistent with temperature variations of a blackbody of ~ 2.5 
X lO^'K that emits in the UV. More recently Pereyra et al. ( I2006I I 
have also found that most of the optical-UV variability observed 
in quasars is due to changes in the UV continuum level which 
can be directly attributed to processes involving changes in the 
disk mass accretion rates. Repeatedly visible spectroscopic ob- 
servations of 2,500 quasars over > 50 day time-frame by Wilhite 
et al. ( 120051 ) have shown that the strongest quasar emission lines 
vary by only 30% as much as the underlying continuum level. 



All these results are therefore consistent with the increasingly 
widely-held belief that the majority of a quasar's UV-optical 
spectrum arises directly from the accretion disk and that changes 
in the level of the underlying UV continua of quasars are proba- 
bly the main driver of the observed variability in these objects. 

4.1. Ensemble Structure Function Dependence 

The ensemble Structure Functions shown in Figure 7 clearly re- 
veal that both the FUV and NUV variability of UVQs are well 
correlated with rest frame time-lags of ^ 300 days, such that the 
UV variability amplitude increases almost linearly with increas- 
ing time-lag. However, for longer time-lags there is a rollover in 
both SF curves, which is also observed at visible wavelengths. 
It is therefore initially tempting to attribute this rollover to a 
change in the underlying physical process that cause the ob- 
served UV variability. However, a rollover in the SF at large 
time-lags has recently been fully explained by the 'damped ran- 
dom walk' model of quasar variability presented by Kelly et al. 
( 120091 ) and MacLeod et al. (120101 ). These authors model quasar 
light curves as a continuous time stochastic process, which re- 
sults in the estimation of a characteristic timescale and amplitude 
for visible flux variations. In the visible, these flux variations are 
thought to result from thermal fluctuations that are driven by an 
underlying stochastic process such as a turbulent magnetic field 
associated with the black hole accretion disk. The shape of their 
power spectrum (i.e. SF) is predicted to be 1 / f^, which then flat- 
tens to white noise at timescales long compared to the relaxation 
time. This relaxation time is interpreted (for quasar flux varia- 
tions) as the time required for the time series of fluctuations to 
become roughly uncorrected. From our present NUV and FUV 
data this relaxation time is ^ 400 days. 

The reader should be reminded, however, that our GALEX 
observations are biased towards the detection of brighter and 
lower-redshift objects which could potentially be making a 
larger contribution to the ensemble SF at long time lags com- 
pared with the higher redshift (and fainter) UVQs. Hence, al- 
though the ensemble SFs shown in Figure 7 are important in- 
dicators of quasar variability, only by discussing the potential 
dependencies shown Figures 8-13 can we better interpret these 
data. 

4.2. Time Lag Dependence 

The NUV and FUV SFs shown in Figures 8 and 9 as a function 
of time-lag follow the same general trend as their ensemble SF 
counterparts shown in Figure 7, in that the amplitude of UV vari- 
ability increases (linearly) as a function of rest frame time-lag 
up to ^ 200 days and thereafter there is a pronounced rollover 
in the NUV and FUV SFs for all redshift values. This overall 
trend is similar, irrespective of the value of UVQ redshift. For 
short time-lags < 20 days there is a definite trend for higher-z 
UVQs to be more NUV and FUV variable than their low- 2 coun- 
terparts, but this trend all but disappears for time-lags in excess 
of 100 days. If we assume that the observed variability is mainly 
due to changes in UV line emission, then we might expect the 
NUV band to exhibit greater variability than that recorded in 
the FUV band for objects with 2 = 0.46 to 1.33 (i.e. the range in 
which the strong Ly-a line enters the NUV band pass). However, 
our observations do not strongly support this notion since the SF 
variability levels of the FUV data for these objects are similar to 
those of their NUV variability. Thus, it seems likely that the ob- 
served FUV variability for the UVQs in this redshift range may 
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be better explained through (stochastic) changes in the underly- 
ing UV continuum level. 

4.3. Luminosity Dependence 

Although it is generally agreed that the intrinsic luminosity of 
a quasar may be linked to matter accreting from a disk onto a 
black hole, the physical mechanism(s) responsible for the ob- 
served changes in variability are still widely debated. As pointed 
out in the Introduction, possible variability mechanisms from 
extrinsic gravitational sight-line microlensing to intrinsic ac- 
cretion disk instabilities have been forwarded to explain the 
observations. Using /. U.E. satellite observations, Paltani & 
Courvoisier ( |1997| l found that quasar variability and luminos- 
ity (absolute magnitude) were related by a power law slope of 
-0.08, such that relative variability decreases as the quasar lumi- 
nosity increases. This anticorrelation between quasar variabil- 
ity and luminosity has also been confirmed at optical wave- 
lengths using large homogeneous data sets (Bauer et al. 120091 
Vanden Berk et al. 2004). However, we note that Ai et al. ( 1201 Oi l 
claim that the inverse dependence of variability on luminosity 
is not significant for a sample of SDSS selected AGN (mostly 
quasars). They claim that their analysis was capable of recover- 
ing AGN with a small variability in the low luminosity regime 
which were largely missed in most of the previous studies. 
Furthermore, previous relationships between quasar luminosity 
and other parameters such as black hole mass and spectral line- 
widths have mostly been derived from single-epoch estimates 
and using (mostly untested) extrapolations from low to high red- 
shift objects (Vestergaard & Peterson ,20061 1 . 

Our present findings shown in Figure 10 indicate that there is 
a general decrease in both NUV and FUV variability as a func- 
tion of increasing absolute magnitude for time-lags > lOOd. The 
opposite effect is observed for time-lags < lOOd (apart from 
the weak exception in the NUV for time-lags < lOd). Figure 
1 1 shows that this general decrease in UV variability as a func- 
tion of luminosity is not strongly coupled to the UVQ redshift, in 
conflict with other (visible) studies that have claimed a redshift- 
variability relationship (Trevese & Vagnetti 2002). Both of the 
NUV and FUV SFs shown in Figure 10 reveal a similar behav- 
ior to those shown in Figures 7 -9 (which were derived indepen- 
dently of luminosity), in that variability amplitudes are greater at 
longer time-lags than for shorter time-lags for all values of QSO 
luminosity. Therefore, it would seem from our present data that 
the intrinsic luminosity of a UVQ may only be weakly connected 
to the origin of the physical mechanism causing the UV variabil- 
ity, and that this mechanism operates over most time-lags and 
is not strongly dependent of the UVQ redshift. As discussed in 
Section 4.1, although the actual underlying physical causes have 
yet to be proven, the application of a 'damped random walk' 
model of stochastic thermal fluctuations from an accretion disk 
can statistically explain the general shape of the ensemble vari- 
ability of quasars, and in particular the shape of their SFs as a 
function of time-lag (Kelly et al. |2009l MacLeod et al. |20T0] |. 
Therefore, in order to explain the majority of our UV variabiUty 
data we do not have to invoke various physical mechanisms that 
operate over different time-lags that depend on various physical 
parameters, but instead the gross behavior of the UV variability 
of quasars can essentially be described by a 1 / f^ dependence 
that flattens to a near constant at low frequencies (i.e. long time- 
lags). This seems to be true apart from the effect noted in Section 
4.2 (and shown in Figures 8 and 9) in which higher-z UVQs ap- 
pear to be more variable than their low-z counterparts for time- 
lags < 20 days. Whatever physical (or observational) effect is re- 



sponsible for this behavior essentially disappears from our data 
for time-lags > 100 days. 

4.4. Radio Loudness - UV variability 

Although the data presented in both Figures 12 and 13 have been 
constructed with relatively few UVQs, we note the near similar- 
ity between the shape of the NUV SFs for both radio loud and 
radio quiet objects. In addition, we note that both Figures also 
show that radio loud UVQs are marginally more variable than 
radio-quiet sources (a factor 2 in the NUV and a factor 1 - 3 in 
the FUV, depending on time-lag), in agreement with the work 
of both Eggers et al. (2000) and Helfand et aU UMTi . Similarly, 
several other visible variability studies (Vanden Berk et al. 2004] 
Ramirez et al. 2009) have revealed a correlation between radio 
loudness and UV/optical variability amplitude that was sugges- 
tive but not conclusive. The source of quasar radio emission is 
thought to originate either in an accretion disk around the central 
black hole or in shock waves in a relativistic jet ejected from a 
region around the central black hole. Since the processes asso- 
ciated with a powerful nuclear radio source probably influence 
UV variability over time-lags 100 days, we might expect the 
greatest difference in variability between radio loud and radio 
quiet objects to appear in this time-regime. Both of our FUV 
and NUV data do not support such a scenario, and to the con- 
trary the (sparsely sampled) FUV SF suggests that the greatest 
variability difference occurs for time lags > 100 days for radio 
sources. Thus, although our data do reveal marginal differences 
between the behavior of radio loud and radio quiet sources, we 
believe that the radio properties of UVQs appear only play a mi- 
nor role in the gross levels of observed UV variability. 



5. Conclusion 

Using archival observations recorded at NUV and FUV wave- 
lengths over a 5+ year timeframe with the NASA Galaxy 
Evolution Explorer ( GALEX) satellite, we have performed an 
ultraviolet variability analysis of quasars that have optical coun- 
terparts in the Sloan Digital Sky Survey (SDSS) DR5 spectro- 
scopic catalog of Schneider et al. (2007). The application of a 
set of (conservative) selection criteria has resulted in the iden- 
tification of 550 quasars in the NUV and 371 in the FUV as 
being statistically significantly variable using data recorded in 
5 or more GALEX observations spread over time periods > 3 
months. Typically we have observed a wide variety of UV vari- 
ability levels ranging from 0.1 1 to 3.0 AB magnitudes arising in 
quasars with redshifts of < z < 2.5. 

In order to explore possible correlations within these vari- 
ability data we have used a Structure Function analysis that 
closely follows the work of Vanden Berk et al. '2004! who 
performed a similar analysis on visible observations of SDSS 
quasars. By subdividing the variable quasars into time-lag and/or 
redshift bins we have been able to obtain a set of (quasi- 
independent) variability relations with respect to the Structure 
Function (SF) for values of rest frame time-lag, luminosity and 
radio loudness. Since our selection criteria for UV variability 
was quite conservative, our analysis of the SFs may only be ap- 
plicable to the more UV variable QSOs. 

The SF - time-lag plots clearly reveal that the amplitudes of 
variability in both the NUV and FUV are far larger than those ob- 
served in quasars at visible wavelengths. Additionally, the levels 
of FUV variability are greater than those observed in the NUV at 
any given value of time-lag. We also find that the amplitudes of 
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both the NUV and FUV variability of UVQs increase as a func- 
tion of rest frame time-lag, irrespective of the value of quasar 
redshift, for time-lags < 200 days. For time-lags > 300 days 
there is a pronounced rollover in the NUV SF for all redshift 
values, which is also observed (with a lower significance) in the 
FUV variability data. These data also show that for time-lags < 
20 days there is an apparent trend for higher-z UVQs to be more 
NUV and FUV variable than their low-2 counterparts. However, 
none of of other plots of SF dependence show any strong de- 
pendence on redshift, and this effect may be due to some other 
physical dependency. Overall, these results are best interpreted 
under the assumption that the UV variability is linked to varia- 
tions in the level of continuum, and not line, emission. 

Our data also indicate that for time-lags > 100 days the 
more luminous UVQs tend to be less NUV and FUV variable, 
whereas the opposite effect is mostly found for shorter time-lags 
(with the weak exception in the NUV for lags of < 10 days). 
This difference in the level of SF variability between long and 
short time-lags is also seen in the ensemble SF - time lag plots, 
which were derived independently of luminosity. We also find 
that UVQ redshift does not seem to have a profound effect on 
any of the variability-luminosity relations (aUitough we note the 
previously mentioned trend for higher- z UVQs to be more NUV 
and FUV variable than their low- 2 counterparts as a function 
of short time-lags). The apparent weak effect of redshift on the 
UV variability-luminosity is perhaps not surprising since lumi- 
nosity is an intrinsic physical property whereas redshift is not. 
These data can possibly be interpreted through the application 
of a 'damped random walk' model of stochastic thermal fluctu- 
ations from an accretion disk which can statistically explain the 
general shape of our VUV SFs ((Kelly et al. 2009, MacLeod et 
al-EOlO). 

Although our data set is small, it seems that radio loud 
sources are marginally more NUV variable as a function of time- 
lag than radio quiet UVQs. Therefore it would seem that what- 
ever the main mechanism is responsible for the observed NUV 
variability in quasars, it seems to operate irrespective of their 
radio flux. Such behavior can also be explained by the damped 
random walk model. Our data also hints that for time-lags > 100 
days radio loud sources may be more FUV variable than radio 
quiet sources, although this conclusion is based on small number 
statistics. 

In summation, our analysis favors a quasar model in which 
UV variability is mainly due to (stochastic) changes in the un- 
derlying continuum level, rather than models that favor gravita- 
tional microlensing or discrete-event processes. Finally, we note 
that as yet, no detailed quasar model cuiTently exists that can 
be directly compared with our observations of NUV and FUV 
variability. 
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Fig. 1. Plot of the number of GALEX repeated visits for quasars 
with DIS integration times of > 1000s and bright quasars with 
AIS observation of > 100s (see text for selection criteria). Full ^ . , . 

hne is for quasars with GALEX NU V observations and dashed ^8' ^' '^.^^'^f '^.'^^^tijjution of Ae UV selected quasars with > 5 

Une is for quasars with FUV observations. u uZT^ lu^ll^'l u^, "'"'^^^^^f "i""^'^'' 

both NU V and FUV observations, while open circles are quasars 

with NUV only observations by GALEX. 
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Fig. 2. Plot of the number UV selected quasars and the time pe- 
riod spanning all of their GALEX observations. Full line is for 
quasars observed in the GALEX NUV pass-band, dashed line 
is for quasars observed in the FUV pass-band. 
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Fig. 4. Plot of the standard deviation (tr) from the weighted mean 

(< m >) of a series of observed NUV magnitudes (m) for the 
sample of quasars with >5 observational visits by GALEX. 
Superposed on this plot (full line) is the theoretical 3-a variance 
from the measured mean NUV magnitude. Only points above 
this variance were selected as being statistically variable in the 
NUV. 



tig. 6. NUV (filled circles) and FUV (open circles) 
3ight curves for the UV-variable quasars (UVQs) SDSS 
1H63915.80+412833.6 (left) and SDSS J221629.33+002340 

(right) repeatedly observed over several years with the GALEX 
satellite. Error bars represent the uncertainty in the measured 
AB magnitude (see text). 
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Fig. 5. Plot of the standard deviation (cr) from the weighted mean 
(< m >) of a series of observed FUV magnitudes (m) for the 
sample of quasars with >5 observational visits by GALEX. 
Superposed on this plot (full line) is the theoretical 3-(t variance 
from the measured mean FUV magnitude. Only points above 
this variance were selected as being statistically variable in the 
FUV. 



Fig. 7. Log-Log plots of the NUV (filled circles) and FUV (open 
circles) structure functions with respect to the rest frame time 
lags for the UV-variable quasars. The data averages are shown 
in 10 bins of equal logarithmically spaced time lags. Best fit 
straight lines (dashed) and quadratic fits (full line) are shown 
for both Structure Functions. Both sets of data reveal a flattening 
of their structure functions for time lags > 300 days. 
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Fig. 8. Log-Log plots of the NUV structure function (separated 
into 3 redshift bins) with respect to the rest frame time lags for 
the UV variable quasars. The data are shown for the same 10 
bins of time lags as Figure 7, but with the data averages plotted 
for redshift bins of size < z < 0.47 (open circles), 0.47 < z < 
1.33 (filled circles) and z > 1.33 (open triangles). The vertical 
lines represent error bars for the binned data points. Linear best- 
fits are also shown for each of the redshift bin data sets; long 
dashed line (open circle data), full line (filled circle data) and 
short dashed line (open triangle data). None of these linear fits 
match the NUV data for time-lags > 300 days. 
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Fig. 9. Log-Log plots of the FU V structure function (in 3 redshift 
bins) with respect to the rest frame time lags for the variable 
quasars. The data are shown for the same 10 bins of time lags as 
Figure 7, with the data averages plotted for redshift bins of size 
< z < 0.47 (open circles), 0.47 < z < 1.33 (filled circles) and 
z > 1.33 (open triangles) . The vertical lines represent error bars 
for the binned data points. Linear best-fits are also shown for 
each of the redshift bin data sets; long dashed line (open circle 
data), full line (filled circle data) and short dashed line (open 
triangle data). None of these linear fits match the FUV data for 
time-lags > 300 days. 
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Fig. 10. Plots of the NUV (left) and FUV (right) structure func- 
tions as a function of quasar luminosity (absolute i magnitude) 
for our sample of UV variable objects. The data are plotted in 
bins of 2 mag width and are shown for four rest frame time-lags 
(t) of Id < t < lOd (crosses), lOd < t < lOOd (asterisks), lOOd 
< t < 500d (open circles) and > 500d (filled circles) for both the 
NUV and FUV observations. Best-fit lines to the data points are 
also shown, revealing that NUV and FUV variabihty amplitudes 
are greater at longer time-lags than for shorter time-lags for all 
values of QSO luminosity 
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Fig. 12. Log-Log plot of the NUV structure function as a func- 
tion of rest frame time-lag for radio bright UVQs (R* > 1, filled 
circles) and radio quiet (R* < 1, open circles) UVQs, as detected 
by the VLA 20cm FIRST survey. Best fit straight lines are shown 
for each of the two data sets. The radio bright NUV data are best 
fit with a Une of slope 0.25 (intercept 0.05) and the radio faint 
sources with a similar slope of 0.24 (intercept 0.08). 
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Fig. 11. Log-Log plots of the NUV (left) and FUV (right) struc- 
ture functions as a function of quasar luminosity (absolute i 
magnitude) for our sample of UV variable objects for rest frame 
time-lags of 1 d < lOd and lOOd < 500d for QSOs with redshifts 
of 0.5 < z < 0.9 (open triangles) and z > 2.03 (filled circles) for 
the NUV data, and redshifts of 0.5 < z < 0.9 (open triangles) 
and lA < z < 2.03 (asterisks). The data points have been con- 
nected by a continuous line for the 1 d < lOd data, and by a 
dashed Une for the lOOd < 500d data for ease of viewing. 
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Fig. 13. Log-Log plot of the FUV structure function as a func- 
tion of rest frame time-lag for radio bright UVQs (R* > 1, filled 
circles) and radio quiet UVQs (R* < 1, open circles) as detected 
by the VLA 20cm FIRST survey. Best fit straight lines are shown 
for each of the two data sets. The radio bright FUV data are best 
fit with a Une of slope 0.16 (intercept 0.22) and the radio faint 
sources with a slope of 0.059 (intercept 0.19). 
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Fig. 14. UV-visible spectrum of a zero redshift quasar showing 
the most prominent emission lines. Also shown are the GALEX 
FUV and NUV passbands together with the redshift (z) values at 
which Ly-a emission enters and leaves each of the passbands 
(adapted from Vanden Berk et al. 2001) 



